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Abstract 

The dynamic scaling of curved interfaces presents features that are strikingly different from those 
of the planar ones. Spherical surfaces above one dimension are flat because the noise is irrelevant in 
such cases. Kinetic roughening is thus a one-dimensional phenomenon characterized by a marginal 
logarithmic amplitude of the fluctuations. Models characterized by a planar dynamical exponent 
z > 1, which include the most common stochastic growth equations, suffer a loss of correlation 
along the interface, and their dynamics reduce to that of the radial random deposition model in the 
long time limit. The consequences in several applications are discussed, and we conclude that it is 
necessary to reexamine some experimental results in which standard scaling analysis was applied. 
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In Nature, one finds a large number of phenomena driven by growing interfaces. To cite a 
few, we can recall wetting fronts, flames, bacterial growth, fluid flow in porous media, thin- 
film deposition.... All these, apparently independent, processes have been studied under 
the same theoretical framework: scaling analysis [1]. Scaling analysis characterizes growing 
surfaces with sets of critical exponents, encoding information about the morphology and 
dynamics of the interfaces. Different phenomena sharing the same critical exponents are 
said to belong to the same universality class. Thus, classifying in terms of universality 
classes allows us to identify how growth dynamics with a very diverse origin can have the 
main physical mechanism in common. Despite its successes, there are some limitations 
associated with scaling analysis so far. For instance, the interfaces are supposed to be 
describable from a Euclidean reference frame, and they do not alter their size during the 
evolution. While these assumptions are reasonable in many cases, on can find as many 
where they do not hold. Biology, for instance, is remarkab 
the approximate spherical symmetry: bacterial colonies 



e for its strong tendency towards 
2], fungi [3], plant calli J], and 



tumors 5j develop rough interfaces that escape the hypothesis of traditional scaling analysis. 
Not only biological systems, but also physical settings such as fluid flow in porous media 1], 
grain-grain displacement in Hele-Shaw cells {(j], adatom and vacancy islands on surfaces 7], 
and atomic ledges bordering crystalline facets present interfaces that either become 

larger as time evolves or have a non-Euclidean geometry. For this reason, this problem has 
been considered as one of the main open questions of scaling analysis [lol ]. 

Traditionally, systems like flames or wetting fronts have been studied in strip geometry 
in well controlled experiments [l| although they appear more commonly as propagating in 
all directions from a small, pointlike region. Even bacterial colonies have been grown in 
this fashion [2J, without knowing the relation between the scaling behavior in this geometry 
and the radial one that appears naturally. Tumors have been an exception to this rule 5| 
since they were grown from a seed and they acquired a radial interface during the evolution. 
However, standard scaling analysis was applied to this case as well, a fact that is being 



increasingly questioned in the literature 
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131 ]. These last two systems are 



particularly important in applications, due to the possibility of understanding and even 
controlling bacterial infections and malignant tumor proliferation. Indeed, a strategy to stop 
tumor growth, based on the planar scaling analysis in j^], was proposed in [ijj]. Following 
this reference, an enhancement of the immune response would be able to slow down tumor 
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growth, and it would eventually lead to the pinning of the tumor interface. However, these 



results enter in contradiction with others present in the medical literature. In [15(, it is 
proposed that the infiltration of a tumor by neutrophils, the most common immune cell, can 
promote tumor progression and the appearance of metastases. Furthermore, in this article, a 
possible therapy is proposed, consisting in inhibiting the chemical interleukin-8, responsible 
for the attraction of the neutrophils to the tumor; this is exactly the opposite strategy to 



that proposed in 141. 

n 

The theoretical study of non-equilibrium radial growth dates back to Eden [16], but the 
use of stochastic growth equations appeared only more recently. To our knowledge, the first 



work in this direction is 



171 ]. where the authors proposed the Kardar-Parisi- Zhang (KPZ) 



equation in reparametrization invariant form. After that, several works appeared studying 



the radial (2D) KPZ equation 18|, |l9j, and the radial (2D) and spherical (3D) Mullins- 



Herring (MH) equation 20|, |2l) . The equation of growth of a general Riemannian surface 
reads 

d t f(s, t) = n(s, t)r[f(s, £)] + $(s, t), (1) 

where the d+l dimensional surface vector r(s, t) = {r Q (s, t)}„=\ runs over the surface 
as s = {s l }f =1 varies in a parameter space. In this equation n stands for the unitary 
vector normal to the surface at r, V contains a deterministic growth mechanism that causes 
growth along the normal n to the surface, and $ is a random force acting on the surface. 
The stochastic term comes from the force in the same way as its deterministic counterpart 
$ = 77.77. The noise r\ is assumed to be a Gaussian variable with zero mean and correlation 
given by 

(77(8, t)r,(s', t')) = ra„(s, t)Ms, t^g-^is - s')S(t - t'), (2) 

where the Einstein summation convention has been adopted. We will now focus on the radial 
geometry r = (r(6, t) cos(#), r(6, t) sin(#)), and will start analyzing the simplest process: 
random deposition. In a Euclidean geometry, the corresponding stochastic growth equation 
is dt h = F + r/(x, t), where F is the constant deposition rate, and rj is a zero mean Gaussian 
noise with correlation given by (r/(x, t)ry(x', t')) = e<5(x — x.')5(t — t'). In this case, the only 
well defined exponent is the growth one (3 = 1/2, which is superuniversal, i. e., independent 
of the spatial dimension. This value can be straightforwardly obtained from the properties 
of white noise: after time integration, we obtain a temporal Wiener process in every spatial 
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point; as its variance increases linearly in time, we obtain the desired result averaging over 
uncorrelated spatial points, this is, the interface width evolves as 

1/2 



Tl 2 



W(L,t) = ( [h(x,t) - h] 



t 



1/2 



(3) 



for all times, where the overbar denotes spatial average over the fixed system size L 



(Li 



Ld), and the brackets the average over different rea 



dimensional radial geometry, the growth equation reads 

d t r = F + r~ 1/2 7](9,t) 
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izations. In the case of a one- 
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(4) 



after the usual small gradient expansion, that assumes \dgr\ <C r [18]. The noise correlation 
is given by (rj(9, t)rj(9', t')) = e8(9 — 9')5(t — t'), but due to its multiplicative nature, we need 
to add an interpretation, which is unambiguously Ito. This is because the multiplicative term 
appears as a consequence of the reparametrization invariance of the correlation, E q .(j2|) , while 
the process has zero mean. The interface width can be defined in this case as Q 

1/2 



W r (t) = ([r(9,t)-f] 2 ) 



(5) 



expression that will be evaluated with the help of a small noise expansion 2l| : r(9,t) = 
R(t) + yfep{9,t), where the noise amplitude e is assumed to be small enough. We get 
R(t) = Ft, for t > t , and 

d t p=(Ft)- 1 / 2 a9,t), (6) 

where r/ = y/e£, that is actually an ordinary stochastic differential equation. The corre- 
sponding Fokker-Planck equation for its probability distribution function is 



d t P{p,t) = (47tFt)- l diP(p,t) 



(7) 



that, by changing the temporal variable r = ln(t), is reduced to the diffusion equation. So 
we see that in this case, the variance is proportional to the logarithm of time, and averaging 
over uncorrelated angular points, we get W r (t) ~ [lnft/to)] 1 / 2 for all times, while for short 
times t > t , W r (t) ~ t 1 ! 2 , a behavior reminiscent to that of the planar case. However, 
the long time behavior shows a marked weakening of the fluctuations, and, as we will show, 
the logarithmic dependence is the footprint of critical dimensionality. The two-dimensional 
random deposition equation is [2l(] 



d t r = F+ (r\sm(9) | 1/2 ) r?(#,<M) 



(8) 



where the zero mean Gaussian noise has the correlation (r](9,<f),t)rj(9' ',<(>' ,t')) = 
e5(8 — d')5(<f) — <j)')5(t — t'). The polar coordinate system r = 

(r(9, (f>) sm(9) cos(</>), r(9, </>) sin(0) sin(</>), r(9, (f>) cos(9)) parameterizes the equator of 
the sphere for 9 = tt/2. In this case, the equation for the small perturbation 
jlp{9 = tt/2,0, t) = r{9 = n/2,<f),t) - Ft is d t p = (Ft)-^(7r/2, <f), t), which Fokker- 
Planck equation 

^P(p,t) = (4 7 r)- 1 (Ft)- 2 ^P(p,t), (9) 

is reduced to the diffusion equation by means of changing the temporal variable r = — f -1 . 
Averaging over the uncorrelated angular points on the equator, we get W r (t) ~ t$ 1 — t^ 1 for 
all times, which reveals a fact that is completely new to this geometry: the width saturates 
in the case of the two-dimensional random deposition process, and the surface becomes flat 
after a transient when some residual roughness develops. This result can be easily generalized 
to any spatial dimension d: the Langevin equation for the small stochastic perturbation for 
d > 2 is 

d t p=(Fty d / 2 aM, (io) 

when restricted to the equator, parameterized by the angle 0. The Fokker-Planck equation, 

d t P(p,t) = (4n)- 1 (Ft)- d d 2 p P(p,t), (11) 

is related to the diffusion equation by the change of variable r = —t l ~ d . So we obtain the 
evolution of the interface width W r (t) ~ t] ) ~ d — t l ~ d for d > 2, revealing that the surface 
becomes flat after a transient that is shorter for higher dimensions. The consequence of 
this fact is very important: d — 1 is the critical dimension for the radial random deposition 
process. Above this dimension, the interface becomes flat, and the one- dimensional setting 
shows critical logarithmic fluctuations. Another interesting feature is that the scaling prop- 
erties of the equator, a one-dimensional cut across the ci-dimensional sphere, are those of the 
higher dimensional manifold in which it is embedded. The differences with the planar case 
come from the fact that radial interfaces grow in all directions. So, while we could picture 
the random deposition process as the growth of uncorrelated columns in the planar case, 
the sources of noise are distributed in a growing area in radial processes. This rends the 
fluctuations weaker for higher spatial dimensions. 



Let us now focus on the radial Edwards- Wilkinson (EW) equation Jja, ll9j 



d t r = Dr~ 2 d 2 e r - Dr~ l + F + r~ 1/2 r/(^, t). (12) 



The small noise expansion r = R + yfep provides us with an equation for the deterministic 
contribution R = F — DR^ 1 , that can be solved to yield 

R(t) = DF- 1 { 1 + W [(FhoD- 1 - 1) exp {fHD' 1 - l)] } , (13) 

where W is the Lambert omega function, and we have assumed R(to) = Fto as initial 
condition and t > to- An inspection of this formula reveals that the value F 2 to = D is 
a fixed point of the dynamics [18]. Initial conditions below the threshold, R(to) < D/F, 
shrink till they collapse, while the supercritical behavior is characterized by an accelerated 
growth rate, that approaches the linear law R(t) ~ Ft in the long time limit. Shrinking of 
subcritical solutions produces a collapse in finite time 

T c = DF~ 2 ]n [D(D - F 2 t y l ] - t . (14) 

In this limit t — > T c , the solutions behave as R — > and R — > — oo. This "bistability", 
expansion vs. collapse, can produce transitions from one of these states to the other when 
noise is present into the system, a possibility that was in explored in [l8| by means of the 
backward Kolmogorov equation and simulations of a lattice model. Indeed, the fluctuations 
can drive part of the interface from a general expanding behavior to the shrinking one, as 
observed in plant calli due to dehydration processes [4|. Despite the interest of the possible 
relation among calli shrinking and the radial EW equation, we will not consider here large 
fluctuations able to change the interface dynamics: rare event solutions of the full stochastic 
equation are not analytic in the small parameter, and are thus suppressed by the small noise 
expansion. We will assume henceforth a large enough initial radius, such that rare events 
able to shift the interface expansion are so infrequent that their appearance can be ignored 
during the evolution time of interest. The stochastic perturbation obeys the equation 

d tP = D{Ft)- 2 {dip + p) + (Ft)- 1 ^, t), (15) 

in the long time limit. Expressing its mean value as a Fourier series (p(6>, t)) = J2 n e%ni Pn(t\ 
we get for the Fourier modes 

p n (t) = exp [DF-\l - n 2 ) (t l - r 1 )] p n (t ), (16) 



reflecting that the n = mode is unstable and the n = ±1 modes are marginal 19) . The 



effect of perturbing the zeroth mode can be estimated by means of the exact solution Eq. ffTBl . 



while the marginality of the other two is possibly related to the deviation from the circular 



shape of the growing interface of the lattice model [18]. The rest of the modes is stable, 
but part of the perturbation persists even in the long time limit, a fact that can explain the 
higher complexity of the radial interface when compared to its planar counterpart |18 ]. All 
of these features are new and due to the non-Euclidean geometry of the surface. We can 
also get an equation for the correlation 



C(9,9',t) = j {nM) (Pn(t)Pm(t)) 



n,m=—oo 



d 

dt 



(Pn(t)p m (t)) 



D 

~FH< 



[2-n 2 -m 2 ) ( Pn (t) Pm {t)) + 



2ixFt 



(17) 



(18) 



and solve it to obtain 
(p n (t)p m (t)) = exp 



D 

FH 

&n,—m 

2ttF 



(m 2 + n 2 - 2) 
" D 



(Pn(t )pm(t ))exp 



D 



(2 - m 2 - n 2 ) 



Ei 



FH 



(2 - m - n 2 ) 



-Ei 



D 



'2 - m 2 



n 



(19) 



where Ei denotes the exponential integral. Assuming an uncorrelated initial condition and 
using the asymptotic equivalence Ei(z) ~ In (\z\) when z ~ yields 



(Pn{t) Pm {t)) 



-ln(f), C(9,9',t) 



5(9 - 9') 



(20) 



when t —>■ oo. This shows that the interface becomes uncorrelated for long times, and it 
is exactly the same result that we would obtain from Eq.([6]). So, except for the bistability 
expansion-collapse, the radial EW dynamics reduces to random deposition for long times, 
a fact that can be easily generalized to an arbitrary dimension. Not only EW, but also the 
simpler radial MH dynamics [3], that shows a monotone linear growth behavior, reduces to 
random deposition in the long time limit, as can be seen upon applying the same argument to 
the —on function Q. The physica. origin of this asymptotic behavior is the diffusive 
(sub diffusive) propagation of the correlations of the EW (MH) dynamics, as given by the 
dynamical exponent z = 2 (z = 4) in the planar case, that, in presence of a linear increase 
of the system size, rends the interface uncorrelated. 
The radial KPZ equation 18j, 119J 

d t r = Dr- 2 d 2 e r - Dr' 1 +F[l + 2- 1 r~ 2 (d e r) 2 ] + r - 1/2 r](9, t), (21) 

retains the first order nonlinearity in the derivatives of the radius accounting for lateral 
growth. The planar KPZ dynamical exponent z = 3/2 shows sub-ballistic behavior, which 



suggests that, assuming the validity of the planar description locally for the radial process 
would yield the random deposition dynamics for long times. Let us show this explicitly. 
Considering just this new nonlinearity in the small perturbation approximation 

d tP = F2- 1 (Fty 2 (d e p)\ (22) 

shows that, upon applying an angular derivative to both sides of this equation and shifting 
the time variable r = (Ft)' 1 , the gradient dep(9,r) obeys Burgers equation. So pertur- 
bations are constant along characteristics in inverse time, implying in turn that they will 
be able to correlate, at most, a finite angle proportional to the size of the perturbation 
and to the inverse of the initial radius. This can be read from the exact implicit solu- 
tion pq = pf^ \9 — peF' 1 ^ 1 — tg 1 )], where pf denotes the initial condition. So, for large 
enough radius t ^> 1, small perturbations will be able to correlate only part of the inter- 
face, that becomes shorter for larger initial times. The hydrodynamic limit will be again 
characterized by a radial random deposition behavior, after a perhaps long and complex 
transient when part of the planar KPZ exponents could be temporally recovered in some 
spatio-temporal scales 

In conclusion, we have shown that the dynamic scaling of non-Euclidean interfaces 
presents strikingly new features. Growing surfaces above d = 1 become flat after some 
transient because the noise is irrelevant in such cases. Only one-dimensional interfaces show 
nontrivial stochastic dynamics, but the noise is much weaker than that of the planar case: 
the random deposition interface width increases only as the square root of the logarithm of 
time. Our results also point to the fact that models showing sub-ballistic propagation of the 
correlations, corresponding to the value of the planar dynamical exponent z > 1, reduce to 
the radial random deposition model in the long time limit, implying the loss of correlation 
of the interface. Taking into account that the most common stochastic growth equations 
fulfill z > 1 in planar geometry, we expect that racfial random dep 08 it,on dynamics will be of 
broad importance to understand the physics of growing radial clusters. In [6| , it is shown an 
interesting effect that appears in radial geometry: the value of the critical exponents might 
be different if we choose a static origin as reference or if we use the border center of mass. 
As noted by the authors, this is due to the nontrivial growth of the border center of mass 
fluctuations. According to our present analysis, the interface fluctuations trivialize in the 
case of higher dimensions, suggesting that this effect is genuinely one-dimensional. This fact 
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can be related to the properties of adatom and vacancy islands on surfaces. These islands 
perform Brownian motion in equilibrium, with their diffusion constants being proportional 
to a negative power of their radius, and their interface dynamics being described by the EW 
and MH equations [7] . Out of equilibrium, adatom clusters might grow by deposition of new 
atoms, while vacancy clusters grow due to evaporation processes. In this regime, it would be 
interesting to observe if the growing islands behave as described by the radial EW and MH 
equations, this is, showing an uncorrelated and logarithmically rough interface. Also, the 
behavior of a rounded crystal facet ledge might be of deep interest in order to understand 
KPZ dynamics. In equilibrium, the roughness exponent is a = 1/3 instead of the random 
walk one a = 1/2, as predicted theoretically and afterwards found experimentally 
This reveals that a nonlinearity in the manner of KPZ is accounting for the equilibrium 
fluctuations. In the nonstationary situation, for growing ledges, we predict a dynamic be- 
havior in the hydrodynamic limit that is totally analogous to that of an island. So this 
physical setting constitutes an excellent test that can clarify whether this hydrodynamic 
limit is unphysically large or experimentally accesible. The implications of our results in 

growing tumors 



medicine can have deep consequences too. The planar scaling analysis o 
can lead to misleading results, and any therapeutic strategy based on it 14J should be taken 
with extreme caution. In light of our results, it seems prudent to reconsider other complex 
growth problems to assess the effects of geometry. As we have shown in the basic models 
considered here, the consequences can indeed be profound. 
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